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Abstract 

We consider the problem of an atomic beam propagating quantum mechanically through an 
atom beam splitter. Casting the problem in an adiabatic representation (in the spirit of the 
Born-Oppenheimer approximation in molecular physics) sheds light on explicit effects due to non- 
adiabatic passage of the atoms through the splitter region. We are thus able to probe the fully three 
dimensional structure of the beam splitter, gathering quantitative information about mode-mixing, 
splitting ratios, and reflection and transmission probabilities. 
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I. INTRODUCTION 



Continuing advances in the production and manipulation of atomic Bose-Einstein con- 
densates (EEC's) are tending toward novel applications in interferometry. EEC's can now 
be produced in situ on surfaces 0,0,0], making them ready for loading into "interferometer- 
on-a-chip" micro-structures. Eeing in close proximity to the chip, the atoms are subject to 
control via magnetic fields generated by wires on the chip. Eecause of their coherence and 
greater brightness, Eose- condensed atoms are expected to improve uponprevious accom- 
plishments with thermal atoms, such as neutral atom guiding, 0, 00, , switching [l^ , 
and multi-mode beamsplittin g |llL ll2l|. Studies of propagation of EEC's through waveguide 
structures are also underway \U 

While the EEC is created in its lowest transverse mode in the guiding potential, keeping it 
in this mode as it travels through the chip remains remains a significant technical challenge. 
For example, it appears that inhomogeneities in the gui ding wires produce field fluctuations 
that can break up the condensate wave function [l^, ll^- Additionally, the very act of 
splitting a condensate into two paths implies a transverse pull on the condensate that can 
excite higher modes. Ideally, the condensate propagates sufficiently slowly that, once in 
its lowest mode, it follows adiabatically into the lowest mode of the split condensate. The 
criterion for this to happen, roughly, is that the condensate velocity in the direction of 
motion be less than Lu, where L is a characteristic length scale over which the beam is 
split, and a; is a characteristic frequency of transverse oscillation in the guiding potential. 
Reference has verified this conclusion numerically, in a two-dimensional model that 
varies the transverse potential in time, at a rate equivalent to the passage of the moving 
condensate through a beam splitter. Populating higher modes can reduce fringe contrast, 
thus spoiling the operation of an interferometer. Diffraction has also been pointed out to 
have negative effects on guiding in general [iG^ . 

Moving too slowly through the beam splitter is, however, potentially dangerous because 
of threshold scattering behavior in a varying potential. In one dimension, a wave incident on 
a scattering potential is reflected with unit probability in the limit of zero collision energy 
jirl l . This same kind of "quantum reflection" will be generically present in beam splitters as 
well, where scattering can occur from changes in the transverse potential as the longitudinal 
coordinate varies. Reflection upon entering the beam splitter region can prove devastating 
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for potential applications such as a Sagnac interferometer. 

Both aspects of instability in an atom interferometer can be expressed in terms of quan- 
tum mechanical scattering theory of the atoms from the guiding potential. Specifically, a 
condensate entering a beam splitter in arm a and in transverse mode m possesses a scat- 
tering amplitude Sam,a'm' for exiting in arm a' in mode m'. In this paper we therefore cast 
the general problem of beam splitting in terms of scattering theory. For the time being we 
restrict our attention to the linear scattering problem, and therefore implicitly consider the 
regime of weak inter-atomic interactions. This is suitable, since the basic question we raise 
is the effect of wave mechanical propagation on the atoms. Note that the weakly interacting 
atom limit is achieved with small atom number, in which case number fluctuations may be 
problematic |18|]. Alternatively, this limit is reached at low atom density, which is achieved 
for a BEC that has expanded longitudinally for some time before entering the beam splitter 
region. 

Restricting our attention to the linear Schrodinger equation opens up a host of powerful 
theoretical tools that have been developed in the context of atomic scattering. In the 
present instance, given the dominant role of non-adiabatic effects, the tool of most use is 
the adiabatic representation. This is analogous to the Born-Oppenheimer approximation in 
molecular physics Specifically, we freeze the value of the longitudinal coordinate z and 
solve the remaining 2-dimensional Schrodinger equation in x-y. The resulting z-dependent 
energy spectrum represents a set of potential curves for following the remaining motion in 
z. This general approach has been applied previously to a model situation in which the 
transverse potential is gently contracted or expanded [la, 120]; here we extend it to realistic 
waveguide geometries. 

This representation has obvious appeal for the problem at hand, since in this level of ap- 
proximation it is assumed that the atoms move infinitely slowly through the beam splitter. 
It is, however, an exact representation of scattering theory, and the leftover non-adiabatic 
corrections, arising from finite propagation velocity, can be explicitly incorporated. We will 
see that non-adiabatic effects have a strong influence on beamsplitters based on experi- 
mentally realistic parameters. The effects of excitation of higher transverse modes and of 
reflection from the beam splitter therefore have a fairly simple interpretation in these explicit 
nonadiabticites. In addition, the successive solution of a set of two-dimensional problems in 
transverse coordinates x — y, followed by a coupled-channel calculation in z, is less numer- 
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ically intensive than than determining the full 3-dimensional solution all at once. Indeed, 
this is why adiabatic representations have found widespread use in chemical physics. Larger 
problems, more closely resembling experimental beam splitters, can therefore be handled. 
This paper is organized as follows: In section II we introduce the model, describing how 
the beamsplitter works in general terms and outlining the theoretical methods used in the 
paper, introducing the main ideas about the adiabatic formalism. In section III we present 
the results obtained from our study, with a focus on the description of the theory itself, and 
how its different components relate to the physics of the problem. 

II. MODEL 

The salient characteristics of a two wire atomic beam-splitter can be realized in the 
following way: a guiding potential is generated by the magnetic field due to two parallel 
current carrying wires and an additional bias field perpendicular to them. By suitably 
decreasing the bias field or the distance between the wires, it is possible to decrease the 
separation between two minima, and thus increase the probability for the atoms to tunnel 
from one to another. 



A. The Beam Splitting Potential. 

We start by considering the magnetic field generated by two infinitely long parallel wires 
lying on a substrate, each carrying a current I in the z direction. Defining the plane of the 
substrate as the x — z plane, we let the z axis lie exactly between the wires, and let y axis 



point to the region above the substrate. 



We then proceed with the addition of two bias fields, one in the z direction, Bbz, and 
one in the x direction, B^x- The first of the two is put in place in order to avoid regions of 
exactly zero field, where Majorana transitions would cause arbitrary spin flips, and therefore 
loss of atoms from the guide. The second of the two fields, when added vectorially to the 
field generated by the wires, generates regions of minimum potential in the x — z plane. 
In particular, for B^^ = fiQl/27id , where /xq is the permeability of free space, and d is the 
separation between the wires, there exists a single potential minimum located on the y axis 
a distance yo = d above the wires. 
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Furthermore, for B^x < B^^ two minima are generated on the y axis, one above and one 
below yo, and for Bb^ > B^^, two minima are again generated above the substrate, but this 
time they are displaced symmetrically to the left and the right of the y — z plane. It is 
this latter regime that we use to generate a beam splitter, letting the wires be fixed, and 
changing the transverse bias field B^^ as a function of z from B^ax to i?mm and back, such 
that BjYiax ^ B^ifi > Bfj^. 



The magnetic field produced by such configuration is therefore |l8, 
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and the consequent potential experienced by the atoms is 



V = gpfiBmpl^l 



(1) 



(2) 



where /ib is the Bohr magneton, Qf is the Lande factor, mp is the total angular momentum 
projection quantum number, and the atoms' spin is aligned with the field at every point in 
space. An example of a guiding potential is illustrated in FigHJ 

The adjustable experimental parameters are therefore the current in the wires J, the val- 
ues of the bias fields B^ax , Bmm and B^z, and the distance d between the wires. Throughout 
this work we choose, for concreteness, d = O.l/im, Bmin = 21. 3G, B^ax = 22.5(7, Bf,z = LOG, 
L = 20 fim, and we let B^xiz) be a fourth order polynomial in z, such that it has zero deriva- 
tive at the center { z = 0) and edges {z = ±L) of the beam splitter. Also we will only 
consider cases in which Bhx{z) reaches its minimum value at z = only, avoiding the char- 
acterization of the trivial evolution of the wave function at a constant field. In particular 
we will consider the following form for the variation of the transverse bias field: 



Bbx{z) = Bmim + {Bmin " ^max) [2(a;/L)^ - (x/L)''] (3) 

Varying L will therefore adjust the adiabaticity of the beamsplitter, whose effects we will 
study in section (III B) 




FIG. 1: Constant Energy surface representing the potential in Eq.©. The parameters chosen for 
this plot are d = 0.1/im, Bmin = 21. 3G, Bmax = 22.bG, B^z = l.OG, L = 20/xm. The surface 
contour is drawn at the energy of the lowest mode of the input arm of the beamsplitter, and we 
define "left" and "right" arms for labeling convenience. 

B. Waveguides as a Scattering Problem 

Because we are going to treat the beamsplitter as a scattering problem, we will begin 
by offering a quick review of scattering theory; in particular we will reproduce the basic 
formulation of the adiabatic treatment of the scattering problem. 

Scattering theory is fundamentally based on the superposition principle, which constrains 
us to the solution of the linear Schrodinger equation. This limit is nonetheless justifiable 
in light of the known problems caused by the interaction between atoms, such as the wave 
function recombination instabilities described in Ref. 2^ . 

The separation between the guides at the input and output ports of the beam-splitter 
is sufficiently great that no tunneling is possible between the guides within the time frame 
of the experiment. The problem is thus divided into two separate regions. We will refer to 
the region \z\ < L as the scattering region. This is the inner region containing the active 
part of the beam-splitter, where all the coupling between the modes takes place. In the 
outer region, defined by \z\ > L,the potential has translational symmetry in z. Solutions 
to the Schrodinger equation in the outer region are therefore trivially found to be products 
of transverse modes and longitudinal plane-wave solutions. The problem is thus reduced to 
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finding solutions inside the scattering region, and match them to the solutions outside to 
find solutions to the global problem. Once these solutions are found it is then possible to 
generate the S-matrix for the system. 

Moreover, since we are matching at the boundary of the scattering region the only infor- 
mation we need is the value of the wave function and its derivative at the boundary, and 
nowhere else. In particular we need to compute the quantity 

J_a^ (4) 
^(r) dn ^ ' 

defined as the logarithmic derivative, where S is the boundary of the scattering region, n 
is the outward normal to the surface S, and \E' is the wave function in the inner region. 
Because the wave function must vanish in the limit of large |x| or \y\, the surface S consists, 
for us, of the two planes z = ±L. To find solutions inside the box we used the R-matrix 
method, formulated in the adiabatic representation. A derivation of this method follows. 



C. The Adiabatic representation 



We start by writing the Shrodinger equation 



-—Vly + V{x,y,z) 



(5) 



2m dz"^ 

with V{x, y, z) as defined in Eq.®. If atoms in the guide were moving infinitely slowly , i.e. 
adiabatically, then the wave function would be well represented by the basis set (j)i{x,y; z) 
with eigenvalues Ui{z), defined as solutions to the equation 



(f)i{x,y;z) = U{z)(j)i{x,y;z). 



(6) 



As in the Born-Oppenheimer approximation, the quantities Ui{z) serve as effective potentials 
for the subsequent motion in z. To recover the effect of finite velocity in z, it would be 
appropriate to expand the wave function in terms of the adiabatic basis in the following 
way: 

^{x,y,z) = Y^Fj{z)(j)j{x,y;z), (7) 
j 

where the z dependence of the coefficient Fi{z) is necessary in order to restore the motion 
in the z coordinate. We should note that the above defined basis functions depend only 
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parametrically on and they are normalized in the following way: 



(f)i{x, y; z)(f)j{x, y; z)dxdy = S, 



(8) 



This normalization implies that all transverse functions must vanish as x,y oo, and 
therefore defines the effective boundary of the scattering region as z = ±L. 

Having defined the basis set we proceed to insert Eqn. ((7j) into Eqn. (0), and subsequently 
project the resulting equation onto (f)i{x,y; z), to obtain the set of coupled equations 



2m 



dz' 



d 



+ UAz)FAz) = EFAz) 



(9) 



where we have defined, as conventional, 



Q 



(t)i{.x,y;z) 



d_ 

dz 



<Pj{x,y;z) 



^^j{x,y;z) 



(10) 
(11) 



P and Q are operators of momentum-like and kinetic energy-like quantities, and thus reflect 
the influence of finite propagation velocity in z. Notice that P and Q vanish by construction 
in the outer region. 

We have thus cast the original 3-dimensional problem into a collection of 2-dimensional 
problems to find (f)i{x,y] z) and Ui{z), and a 1-dimensional coupled channel problem to find 
Fi{z). The advantages of this shift in paradigm are twofold: on the one hand, a very com- 
plicated and computationally lengthy problem is turned into a simpler and computationally 
manageable problem; On the other hand, the adiabatic approach lends itself very naturally 
to approximations and qualitative understanding of the underlying physics. 

D. The R-Matrix Method 



As mentioned earlier, solving the scattering problem implies finding the logarithmic 
derivative b as defined in Eq.Q. In atomic structure physics, there is a well known 
variational principle for the logarithmic derivative, which follows rather simply from the 
Schrodinger equation 



(12) 



where denotes an integral over the volume of the scattering region, while V is the scat- 
tering potential, and is a surface integral over the surface bounding the scattering region. 

The typical approach to the problem, at this point is to expand the wave function in 
a complete set of basis functions Y^., to get z) = J2kCkyk{x,y, z), and take matrix 

elements of the operators in Eon. (|12p with respect to such basis, to obtain the following 
generalized eigenvalue problem [23|: 

tC = bkC, (13) 



where 



A,; 



VY:-VY, + ^Y*iV-E)Y, 



(14) 
(15) 



This is the form that the eigenchannel R-matrix takes in a diabatic representation. 

The solutions of the eigenvalue problem consists of a set of eigenvectors C"^, and a set 
of eigenvalues 6*^, representing the logarithmic derivatives of the functions = J^kC^Yj^. 
The newly introduced index a refers to the different possible internal states of the system, 
called R-matrix eigenchannels. The concept of eigenchannel in scattering theory can be 
understood by analogy with the concept of eigenstate in bound state problems. In fact, 
just like the energy variational principle leads to an eigenvalue problem for bound state 
eigenfunction and corresponding energies, the variational principle in Eq.()12|) leads to a set 
of eigenchannels, with corresponding eigen-logarithmic derivatives. 

As we mentioned earlier, we plan to work using the adiabatic basis defined in Eq.Q, so 
we expand Eq. p2|l in terms of this set, and obtain the following variational principle: 



+ '^F:{U,{z) - E)F,} + 2F*P,,§-F, + i^*Q,,F, 



(16) 



Since the adiabatic basis only defines motion in the transverse coordinates, it remains 
to expand the longitudinal functions Fi{z) with an arbitrary set of z dependent functions, 
in our particular case we chose basis-splines, in the form Fi{z) = J2k (^ikUki^) ■ We can now 
write the adiabatic equivalent to Eq. (|T!?jl . 

In order to simplify the notation, we combine the indices i, k into the index a, so that Cq 
becomes a vector, and we write 

t^c = 6A°c, (17) 



where 



Kp = I dz 



d 3 ^TTi I d d 

g^y*c9^yp + -^Vai^iiz) - E)yp \ 6ij + ylPij—yp + -^ylPjiVp + ylQijy0 

= y:{^)yp{mr (18) 

In the above equations we have written the P- matrix portion of Eq. p6p in an Hermitian 
form by integrating by parts, and setting the resulting surface integral to zero, using the 
fact that by definition all couplings P, Q must vanish outside the scattering region. 



E. The Outer Region: Matching and Physical Consideration 

Having solved Eq. ()17|) . one obtains a set of eigenvalues 6°", and a set of eigenvectors c'^. 
It therefore follows that on the boundaries S of the scattering region we can connect the 
inner and outer solutions by: 



vi/'^(x, S) = ^ F^{E)<f),{x, y, S) = ^ c^y,(S)0, (x, y, S) 

j a 

(19) 

where ki = ^2m{E — f/i(S)) is real for E > f/i(S), and imaginary for E < f/i(S), and 
S = ±L. At a particular incoming energy, we define a channel with real ki to be "open" 
(meaning energetically available) an channel with imaginary ki to be "closed" . If a channel 
i is closed we set Aj = 0, to avoid unphysical divergences. A similar argument is valid for 
the derivative of the wave function: 

^vI/-(x,y,S) = -6'^vl>-(a;,2/,S) 

(20) 

Eqs p9|2()j) . together with the orthonormality of the set (pi, and the assumption of unit 
incoming flux, imply that -fJ(S) and its derivatives can be written as a linear combination 
of the form: 
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at; 



< 



(21) 





The quantity is a factor which serves to connect the normahzations of the two equa- 
tions. On the other hand Sij is the scattering matrix of the system (often referred to as 
S-matrix), and it represents the probabihty amphtude to enter the beam sphtter in channel 
j, and exit it in channel i, or vice versa, since 5* is Hermitian due to time reversal symmetry. 

Moreover, since the equation is true on the whole of the boundary, the channel index 
describes the probability amplitude for the atom to be found at either end of the beam 
splitter (in fact at any particular arm of the beam splitter), in some particular mode. This 
allows us to calculate mode mixing, as well as reflection and transmission amplitudes. The 
above system of equations can be solved for the unknowns 5*,^ and Nj . 

F. Solving the Equations: Considerations on Numerical and Mathematical Details 

The numerical problem consists of two main parts. The first is to find the transverse 
eigenmodes (f)i{x,y] z). This is accomplished by solving Eq.® at various values of z in such 
a way that the adiabatic curves Ui{z) may be interpolated easily. We accomplish this task 
by generating a Hamiltonian matrix, again using b-sp lines as a basis set, and diagonalizing 
it at various values of z. 

Furthermore one needs to evaluate the P and Q matrices in Eqs.( fTn^llj) . To do this one 
may exploit the Hellmann-Feynman theorem to obtain the following expressions: j^] 




(22) 



and 



Qii 



E 



ki 



(23) 




where 




(24) 



11 



We adopt a common approximation whereby Qij = for j ^ i. The second part consists 
of a scattering problem on the adiabatic curves, by choosing a basis set yi{z). For our 
calculations we use b-splines. 0, 

The guiding potential V{x, y, z) in Fig. ^ exhibits a reflection symmetry about the x — y 
plane. Such a symmetry implies that there is no coupling between even and odd transverse 
modes of the beam splitter. This in turn implies that by describing the problem in a basis 
of even-odd modes it is possible to solve two smaller problems, significantly reducing the 
computational effort. At the end of the calculation it is then possible to perform a change of 
basis to a "left-right" set describing the "left" and "right" arms of the beamsplitter, where 
"right" = "even-odd", and "left" = "even+odd." 



III. RESULTS 



Having described the general formalism, we proceed to report some quantitative results. 
In particular, we use the parameters described in the caption of Fig. and study the 
behavior of the system as we vary the length L over which the beam is split. We focus 
especially on the non adiabatic characteristics of the beamsplitter, namely reflection and 
higher mode excitation. 

The parameters that generate the guiding potential in our model are consistent with those 
in recent chip-based experiments [l^i ll^ . The major difference is that our model guides 
lie close to the substrate, thus tightly confining the atoms in the transverse direction. At 
reasonable atom velocities of several cm/sec, only two modes are then energetically open, 
simplifying the calculations and interpretation in this pilot study. More realistic beamsplit- 
ters can be handled by including the appropriate number of modes in the calculation. 



A. The Adiabatic Curves 



The simplest level of approximation for the problem is to consider only the first even 
and odd mode of the structure, and, analogously to the Born-Oppenheimer approximation, 
ignore all higher modes and couplings. Within the framework of such an approximation we 
see that the Born-Oppenheimer potential depends only on the transverse frequency of the 
guide, which is highest at the entrance and exit of the beamsplitter and lowest in the center. 
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giving rise to curves resembling smoothed square wells. As it turns out the predictions of 
this simple model prove to be grossly inadequate when compared to full coupled channel 
calculation. The reason for this is that the Born-Oppenheimer channels are strongly coupled 
by nonadiabatic effects. 

To suggest how big a correction nonadiabatic effects are, we compare the lowest-lying 
Born-Oppenheimer potential Uo{R) (dashed line in Fig. 2) to the so-called "adiabatic" 
potential, defined by Uq^^ (z) = Uo{z) + Qoo{z) (solid line). The Qoo{z) term represents an 
effect of the transverse momentum on the longitudinal motion. As the guiding potential 
varies as a function of z, the paths of the atoms follow the centers of the guides. This 
causes the atoms to acquire transverse momentum, which removes kinetic energy from the 
longitudinal motion. Thus Qoo is a positive correction. 

In chemical physics applications, the adiabatic curve is sometimes, but not always, a 
better single-channel representation of the problem 27|. In our case, it usefully incorporates 
a primary effect arising from nonadiabaticity. Namely, Uq'^'^ possesses a barrier at the input 
of the splitter. This barrier reflects the fact that kinetic energy spent in transverse motion 
halts motion in the longitudinal direction. Effects of this barrier are evident in the fully- 
converged scattering calculations, below. 

A more complete set of effective adiabatic curves for the first few even and odd modes is 
shown in FigOl For kinetic energies greater than ~ 20/iK, excited state potentials are ener- 
getically allowed in the scattering region. The corresponding mode mixing can be thought of 
as the "sloshing" of the condensate as it is pulled side to side in the potential. Even if these 
excited channels are not energetically allowed, they may (and do) still perturb propagation 
in the lowest mode. Since the length L of the beamsplitter is, in our model, thousands of 
times larger than the longitudinal de Broglie wavelength of the atoms, even a small coupling 
between channels can cause a drastic change in phase shift. This implies that we need a 
fully coupled channel calculation to solve the problem quantitatively. 

Channel coupling is achieved through off-diagonal elements of the P matrix, several of 
which are shown in Fig|3 As expected, the couplings Poi between the lowest channel and 
higher channels i diminishes as i gets larger. Also, as implied by Eq. ()22|) the coupling is 
strongest where the potential is steepest in the longitudinal direction (i.e. z = ±1 in the 
figure). 



13 



-2-1012 

z((im) 



FIG. 2: Comparison of lowest-lying Born-Oppenheimer (Uq) and adiabatic (Ug ) curves, for the 
beamsplitter in Fig. Q 
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FIG. 3: Effective adiabatic potentials for the configuration in Fig. ^ Each curve corresponds to a 
different transverse mode of the beamsplitter. Because of the intrinsic symmetry of the potential, 
even and odd modes can be treated separately. The first six even and odd modes of the structure 
are depicted here. 

B. General Features of Scattering 

Having defined the terms of the problem, and calculated the adiabatic curves and cou- 
plings, we solve the scattering equations and extract the S-matrices of the system. All figures 
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FIG. 4: P-matrix elements coupling the first transverse even and odd modes in Fig. |21to selected 
higher ones. 

shown to this point refer to a beamsplitter with L = 20fim, which is one in which most of the 
typical features are present. Fig. |31 shows the absolute values of selected S-matrix elements 
for this configuration, which represent the probabilities for various outcomes. In particular 
we show the probabilities to exit in the various arm of the beam splitter, assuming unit 
input flux from the left arm of the splitter, as defined in Fig. ^ At the incident energies 
shown in Fig. only the lowest mode in each arm is energetically accessible. This typical 
case is illustrative of the basic elements of the beamsplitter. 

In this beamsplitter the largest probabilities (dashed and short-dashed lines in Fig. 
correspond to transmission, with the probability alternating between left and right arms. 
Thus approximately 50-50 beamsplitting is possible at atom energies where these two curves 
cross. Moreover, the sum of the left and right transmission probabilities is almost, but not 
quite, equal to unity. This can be seen in the slowly decreasing reflection probabilities (solid 
and dotted lines) in the figure. The general features of beamsplitting are preserved under a 
convolution in energy, as exhibited in Fig. Elb). Here and in what follows, convolution is 
used to simplify the appearance of the calculations. 

The reflection probabilities also exhibit a similar left-right oscillation as a function of 
energy. In addition, they exhibit a much faster oscillation. This faster oscillation is familiar 
from one-dimensional scattering from a potential, with one oscillation being added each time 
the energy increases enough to introduce a new de Broglie wavelength into the scattering 

15 



10 




10 





\ 1 / 


\ ' '/ \ 1' /~\ ">-^-"-'\^-'-- i' ' 




— Reflection Left 






Reflection Right 






— Transmission Left 






— Transmission Right 






1 , 1 


I.I.I 




4 6 


8 10 12 



V (cm/s) 



FIG. 5: Probability for exit in the various arms of the beamspHtter depicted in Fig. ^ versus 
incident atom velocity. This figure assumes that the atoms have entered in the lowest mode of 
the left arm. The velocities shown correspond to atom energies below the second threshold, thus 
suppressing higher mode excitations. The bottom panel is a 26fj,K Gaussian convolution of the 
one above. 

region 2^. Here the oscillations are numerous, since the guiding potential is thousands of de 
Broglie wavelengths long. (These oscillations are of course also present in the transmission 
probabilities, but are too small to be seen on the scale of the figure.) 

For smaller values of L , the beamsplitter is badly non-adiabatic, and even qualitative 
features of beamsplitting fade. Fig. IHl shows such a non-adiabatic case, with L = l/im. 
The effect of the input barrier described in Fig. |21 is now much larger, suppressing all 
transmission up to input velocities of about 5cm/s. As the kinetic energy reaches the energy 
of the barrier, the probability exhibits resonant behavior by the presence of spikes in the S- 
matrices. Though mostly washed out by convolution, these features would in principle cause 
transparency of the barrier at extremely well defined velocities, where the kinetic energy 
equals the energy of a metastable boundstate. At higher atom velocities, above the input 
barrier, reflection remains extremely likely, and even the basic action of the beamsplitter is 
destroyed. 
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FIG. 6: Scattering probabilities as in Fig. [3 but for the extremely non adiabatic case L = Ifim. 
Only the convolved plot is shown, with a width of 26 /xK. 

C. Towards Adiabaticity 

Fig. [7| shows reflection probabiUties versus atom velocity, for various values of the beam- 
splitter lengths L. These results are convolved over an energy width of IGfiK, to emphasize 
the overall probability rather than the oscillatory structure. For L > 2/im, reflection de- 
creases nearly linearly on this semi-log plot, suggesting an exponential decrease of reflection 
probability with velocity. Reflection also decreases with increasing L, as expected for an 
increasingly adiabatic beamsplitter. The features noticeable around 12.5cm/s and 13.5cm/s 
represent cusps at the thresholds for the second and third mode to become energetically 
available, smoothed out by convolution. 

Finally, in Fig. IHlwe plot the total transmission to modes higher than the flrst, for input 
velocities higher than the second mode threshold. As might be expected, the probability 
to generate higher modes grows as a function of atom velocity. Countering this trend, the 
probability again diminishes as the length L becomes longer. 

IV. CONCLUSIONS 

We have developed a novel approach to the analysis of non-interacting atomic beams 
traveling through waveguides, based on the adiabatic representation of scattering theory. 
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FIG. 7: Total reflection probabilities for multi-mode beamsplitter of different lengths. Prom top to 
bottom 1/ = 1, 2, 5, 7, 15, 30^m. The cusps around 12.5 and 13 cm/s are the effects of the second 
and third thresholds becoming energetically available. The plot represents a IQfiK width Gaussian 
convolution. 




V (cm/s) 



FIG. 8: Total transmission to higher modes for different lengths beamsplitter. From top to bottom 
L = 1, 2, 5, 7, 15, 30/xm. The plot represents a 3.25f^K width Gaussian convolution. 

This method, originally developed for the study of molecular collision theory, is known to 
be very flexible, and could be applied to many other guiding geometries. We applied this 
approach to the study of a two wire atomic beam-splitter, both to illustrate the method and 
to explore a particular guiding geometry. We have found that the nonadiabatic couplings 
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play a significant role. Because we have deliberately focused on a tightly-confining geometry, 
it is likely that nonadiabaic effects are even more significant in realistic beamsplitters. This 
will be a topic of future study. 
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